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Abstract. In this paper, we use some Fourier analysis techniques to find an 
exact solution to the Cauchy problem for the n-dimensional biwave equation 
in the upper half-space K" X [0, +00). 



1. Introduction 

The Cauchy initial value problem for the n-dimcnsional biwave equation consists 
in finding a scalar function u G C 4 (R n x [0, +00)) such that for (x, i)eK"x (0, +00) 
then 

- a 2 (JL - 6 2 a) u(x, t) = f{x, t), a 2 >b 2 > 0, (1.1) 
together with the initial conditions 

u(x, 0) = fa(x), ^g(x, 0) = <Mx), — (x, 0) = -^(x, 0) = fa(x), (1.2) 

for [x,t) ePx {0}. 

The biwave equation has been studied in some models related to the mathe- 
matical theory of elasticity. Let us consider the mathematical formulation for the 
displacement equation of a homogeneous isotropic elastic body. Remark that, the 
Newton's second law leads to the Cauchy's motion equation of an elastic body, 
which takes the form 

V-<r + f = pu, (1.3) 
where ex is the Cauchy stress tensor field, u is the displacement vector field, f is 
the vector field of body force per unit volume and p is the mass density. 
The infinitesimal strain tensor field is given by the equation 

e= | [Vu + (Vu) T ] . (1.4) 

Moreover, the Hooke's law for homogeneous isotropic bodies has the form 

er = A tracc(e) I + 2/i £, (1-5) 

where A, /i > are Lame's parameters and I is the second-order identity tensor. 
Substituting the strain-displacement equation (|1.4j) and the Hooke's equation (|1.5|l 
into the equilibrium equation (|1.3[) . we obtain the Navier's elastodynamic wave 
equation 

(A + ^i)Vdiv(u) + /iAu + f = p ii. (1.6) 
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This equation in the Cartesian coordinates has the form 

, d I du j \ . „ d 2 Uk , ■, — 




Let us denote a 2 = (A + 2/i)/p, 6 2 = p/p, then (|1 .€>[) can be rewritten as 
d 2 \ 

— -b 2 Aju- (a 2 - 6 2 )Vdiv(u) - f/p = 0. (1.7) 

It is easy to show that the equation (|1.7p has a solution in the following form 

d 2 



x , a 2 A)w+(a 2 - 6 2 )Vdiv(w), (1.8) 
at 



where w is a solution to the biwave equation 



f) 2 \ / B 2 

° .2 a i / l2. 



dt 2 aA ){w- bA )- = i/p - (L9) 

This formula is called as Cauchy-Kovalevski-Somigliana solution to the elastody- 
namic wave equation. Indeed, substituting (|1.8[) - ()1.9[) to the left-hand side of (|1.7p . 
we have 

C = ( * - b 2 A) ((*- a 2 A) w + (a 2 6 2 )Vdiv(w) 



dt 2 J \\dt 2 
dt 2 



-{a 2 - 6 2 )V ( ( |^ - a 2 A I div(w) + (a 2 - fe 2 )Adiv(w) I - f/p. 



Note that 



- a 2 Aj div(w) + (a 2 - 6 2 )Adiv(w) = ( - b 2 Aj div(w). 
Therefore, we get 

b2 A - w + (« 2 - b2 ) - &2a ") vdiv(w) 



\dt 2 J \dt 2 J x ' \dt 2 

-(« 2 - &2 ) V (^5 - 52 A ) div ( w ) - { /P = °- 

For more explanations about physical context, we refer the reader to (HI04, 
IMuslOl ISom50j . Actually, there are not many mathematical papers related to bi- 
wave equations because it gets more difficult when studying high-order PDEs. In 
some recent researches, the symmetry analysis of biwave equations is considered 
and exact solutions are obtained by Fushchych, Roman and Zhdanov FRZ 96] ; the 
existence and uniqueness of the solution to Cauchy initial value problem, bounded 
valued problem are given by Korzyuk, Cheb and Konopelko [KKClOl IKC07] ; the 
finite element methods for approximations of biwave equation are developed by 
Feng and Neilan FNIO] IFNllj . In our present work, the main result is to show the 
exact classical solution to the Cauchy initial value problem for the n-dimensional 
biwave equation by using some techniques of Fourier analysis. 

Returning to the Cauchy problem for the biwave equation (jl.ip , we suppose that 
^0:01:02,03, and / are elements in Schwartz space <S(R n ) of rapidly decreasing 
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functions on K n . Remark that, an indefinitely difierentiable function <f> is called 
rapidly decreasing when (j> and all its derivatives are required to satisfy that 



a.fi = SUp 



d_ 

dx 



4>(x) 



< oo, 



for every multi-index a and /3. The Fourier transform of Schwartz function <fi G 
5(m n ) is defined by 

F[4>] (0 = ^(0= / e- ix -t<j>(x)dx. 



The convolution of two integrable functions (j) and if) is written as It is defined 

as the integral of the product of the two functions after one is reversed and shifted. 
As such, it is a particular kind of integral transform: 



V0(*)= / <t>{r)^{t-r)dr. 



In the Euclidean space M n , the spherical mean of an integrable function cj> around a 
point x is the average of all values of that function on a sphere of radius R centered 
at that point, i.e. it is defined by the formula 



M R (4>)(x) 



1 



fay)da(y) 



1 



(x + Ry)da(y), 



W n R n x JdB(x,R) "' "' w « JdB(0,l) 

where u> n is the surface area of the rt-dimensional unit ball and a is the spherical 



measure area. 



2. Main results 

The Cauchy problem for the homogeneous biwave equation in R" x [0, +00) that 
we will be studying in this section, reads as follows 

d 2 o.\ fd 2 - 



dt 2 



a z A 



8t 2 



b 2 Au 



with the initial conditions 



(hi 



0o (x) , Tt 



01 {x) ' w 



t=0 



0, a 2 >b 2 > 0, 



Or il 

fa (x) , 



(2.1) 



fa (x) , (2.2) 



where 



>0,<Pl 



ifaifa are Schwartz functions. 



The equation (|2.1j) can be rewritten as a fourth-order PDE, which has the following 
form 



Taking Fourier transform to the both sides of the equation (|2.3p , we obtain 

W u(Z,t) + (a 2 + b 2 m 2 Q^u(^t)+a 2 b 2 \^\ 4 u^t) = 0. 
This fourth ODE has the general solution, which takes the form 

u (£, t) = d cos (a\£\t) + C* 2 sin (a\£\t) + C 3 cos {b\£\t) + C± sin {b\£\t) , 



(2.3) 
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where parameters Ci, C2, C3, C4 are determined from the initial conditions: 

du(H,t) 



«(e,*)it= = ^1+^3 = ^0 (o. 



at 



aC 2 \£\ + bC^\ = <\> x (0, 



t=o 



dt 2 



-a 2 Ci|e| 2 -6 2 C 3 |e| 2 =0 2 (0, 



t=o 



d 3 u(£,0) 



dt 3 



Solving above system of equations, we easily get the image of solution u{x, t) via 
Fourier transform given by 

u (£, i) = — r^T77T5 cos (a\Qt) — — sin (a|^|t) + 



(a 2 -6 2 )|£| 2 



(o 3 - a6 2 ) ICI 3 



cos (6|Cj*) + — rn — i-n i^is — sm ( 6 KW 



(a 2 -6 2 )|£| 2 
or by the rewritten form 



(a 2 6-6 3 )|e| 3 



6 2 — a 2 — 



a 2 — 6 2 a 2 — 6 2 

1 (£) FT] + 7T7 [5\« (?) 



^3 (Q 
a 2 -6 2 



a(a 2 - 6 2 ) 



|f| 6(a 2 -6 2 ) 



lsin(6|£|t) 1 sin. (a|£|f) 



6 Id 3 a Id 3 



02 (Q 
a 2 -b 2 



cos(&|£|i) cos(a|£|t) 



kl 2 



(2.4) 

In the next sequence, we will find the inverse formula of (|2.4[) and obtain an exact 
solution to the equation (|2.1[) . 

Theorem 2.1. The Cauchy problem for the homogeneous biwave equation in M x 
[0, +00) has the following solution 



u(x, t) 



2ab(a 2 - b 2 ) 



x-\-at x-\-bt x — bt y 

-b 3 J Mv)dy + a 3 J 0i(y)dy-afe J J cj) 2 (u)dudy+ 

x — at x — bt x — at 



x+at y x+at y r x+bt y r 

+ab J J (fa(u)dudy + b J J J foi^dwdrdy — a J J J foi^duidrdy 

x+bt x-at x-bt 



-a& 3 o (x + at) - ab 3 <j)o (x - at) + a 3 b(j) (x + bt) + a 3 b<j) (x - bt) 



(2.5) 
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Proof. We have some facts that 

e »a|£|t _|_ e — ia|f|t gW|fl* _ g— ia|f|t 

cos(a|£|t) = , sin(a|£|i) = . 

at 

sin _ e w l«l* - e - * ^' _ 1 



lei 2^| 



-at 

at y at y 





at y r 

sin(a|f|i) e M l«l* - e- ia l«l* I f f f M \ UJ , . a* 

TF^ " SEP " "2 J J J e K ' dudTd » + W 

-at 

Moreover, 

— oo 

where 5 (x) is the Dirac delta function. Hence, by the property of the Dirac's delta 
function, we note that 

7 [(6 (x + at) * O (x))] + T [(6 (x - at) * O (x))} 



^ (0 



and 



/at \ at 

| e W^j=I j F[5{x + e)*<t>i{x)]M, 



— at / —at 



at y at y 

<>, U) ( -\ f fe^dudy-^ j fe-^dudy 



o o 

at y at y 



■— J J J- [5 (x + u) * cf>2 (x)} dudy — — J J ' T '[6 (x — u) * 02 {%)] dudy 

J 

(at y t \ at y t 

-\ j j J e^dudrdy) = -1 y j j T[5{x + u)*fo{x)\dudTdy. 
-at / -at 

Substituting the above identities into the formula (|2.4[) . we obtain that 

6 2 J" [<5 (x + at) * O O)] + ? [5 ( x - at ) * 0o {x)\ 



{a 2 - b 2 ) 2 
a 2 T[S(x + bt) * 0o (x)} + J 7 [8(x- bt) * O (x)] 



(a 2 - b 2 ) 2 

at bt 

— at —bt 

(at y at y 

~J J T [5 {x + u) * 02 (x)] dudy - - J J ' I '[6 {x - u) * 2 (x)] dudy 
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bt y bt y 

J- [S (x + u) * <p2 (x)] dudy — — I I J- [6 (x — u) * 02 (x)] dudy 



{a 2 -b 2 )\2jj y 2 





at y t 

1 1 



(a 3 - ab 2 ) 2 
1 1 



J- [S (x + u) * 03 (a;)] dudrdy 



-at 
bt y t 



(a 2 b~b 3 )2 

-bt 



J- [S (x + u) * 03 (x)] dudrdy. 



Consequently, we get the inverse formula of u given by 

b 2 (5 (x + at) * 0o (a:)) + (6(x- at) * O (a:)) 



u (x, t) 



{a 2 - b 2 ) 2 
a 2 {5 (x + bt) * 0o 0*0) + (S(x- bt) * O (x)) 



(a 2 - b 2 ) 2 

at bt 

— at —bt 

(at y at y 

~J J (S (x + u) * 02 (x)) dudy - i J J {5 (x - u) * 2 (x)) dudy 


(bt y bt y 

-i J J (5(x + u)* 02 (x)) dudy - i J J {6 (x - u) * 2 (a;)) dwdy 


at y t 

(S (x + u) * 03 (x)) dudrdy 



1 1 

(a 3 -a6 2 ) 2 

1 1 



-at 
bt y t 



(a 2 6 - 6 3 ) 2 



((5 (x + u) * 03 (x)) dudrdy. 

-bt o o 

The last formula is equivalent to the one given at ()2.5[) , so the theorem is proved. □ 



For the generalized case, we will use the following result: 

Lemma 2.2. For an odd number n > 3, m — and < k < to then 

- ' 2 ~ ~ 

f R e -ISI (jR 2 _ s 2 yn-k ds = 1 ( 1 d\ k ( _1 I e -i*-t da{x) ] . 
J-r 2*kl\RdRj yu^RJsB^ yi ) 

For the proof, we refer the reader to Torchinsky's paper in |Tor09j . Note that, 
in the case k = to, it follows that 



sin H^) = i r e ^\ ds = _l_ _J_ / e -i,.i M 

l£l 2J-at 2™+ 1 ml\a 2 tdt) \uj n -iat Jg B (0,at) 



(2.6) 
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Differentiating with respect to t, then 

cos(a|f|i) = - } , j- f 4-1-1 I — " — / e- ix -«d<j(a:) I . (2.7) 

v 1"=.! > 2 m + 1 m\adt \a?t dt ) \ u} n -iat 1 ' 



dB{Q,at) 



On the other hand, we have 

/• a / 



COS 



(TO co S (a m _ /-slaWel) 1 . , ^ 



' f-l-1 I I e-^daix)) ds. (2.8) 

Vs^sy \U)„-lS J dB( o„) ) 



Ibt 2 m + 1 ml\sdsj yuJn-is J dB ( ,s) 
Integrating the above formula with respect to t, it implies that 
1 sin(6|£|t) 1 sin(a|£|t) 1 



b |C| 3 a ICI 3 2 
1 /l 9\ m / 1 



a v 



e lr ^drdsdiy 

Jbv J-s 



, ,(-£-) l—^—f e- lx ^da{x)\dsdv. (2.9) 
h Jb„ 2 m +im! \sds) \u n - lS J eB{0iS) K 'j 

Moreover, for each function 9 G <S(R"), we also have 

A- |e| V / 2™+im! Va 2 ^t; \u n -iat J dBi0M) J Rn VW W 



,2, 



Wn-lO* JdB(0,at) ■ 

9(x)da(x) 



2 m+1 m\ \a 2 tdtj \u; n -iat JdB(o,at) 
Therefore, we conclude that 



2m+i m \\aHdtJ \u n -iatJ aBi0M) y 'J \£\ 
By the Fourier inversion and convolution formulas, we obtain the identity 

1 f :rv^ sin ( a l£l*) ^.z 



(0 YJr' e*-'d£ 

0i(y)dcr{y) 



(27r)«y R „™ \z\ 
i / 1 d\ m ( i 



2 m + 1 m! \aHdtJ yu^at J 9B{x at) 
i / i \ m 

^M)l(si) (<->-*-(*.)(.)). 

Applying the same way for the expressions (|2 . T[) - (j2 .9[) and substituting the ob- 
tained identities into the formula (|2.4[) . we have found an exact solution to the 
n-dimensional biwave equation, where n > 3 is an odd number: 

Theorem 2.3. The Cauchy initial value problem for the homogeneous n-dimensional 
biwave equation, where n > 3 is an odd number, has the following solution 



, N 1 

u[x, t) 



(n-2)!!(a 2 - 6 2 ) 



b dt V bH dt 
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n—3 

1 d 



-(^1) (•W)"- 2 .U W (o,),.r.)-^(-l 7 ^) ((«/.."- 2 .U„.o l „,)) + 



i a 



i 2 t a< 



s 9s 



((«) n - a M.(^)(a!))d8 



./ bv 



! M s (</> 3 )(a;))ds^ 



Now we consider the case when n is an even number. The Hadamard's method of 
descent (see e.g. |Had53] ) is useful to connect with the case in the odd dimensional 
space R n+1 . For fixed T > 0, we choose a Schwartz function rj £ <S(K), such that 
r)(x n +i) — 1 for all |x„+i| < nT . Let us denote 

4> i (x 1 ,x 2 , ...,x n ,x n+1 ) = 4> i (x 1 ,x 2 , ...,x n )r)(x n+1 ), i = 0, 1,2,3. 



It is easy to see that . 



n+l 



). For \x n +i | < T,t < T, the solution u(xi, X2,-., x n ^ 



to the Cauchy problem for n + 1-dimensional biwave equation with initial valued 
functions <fii, i = 0, 1, 2, 3 does not depend on x n +i- In particular, 

u(xi,X 2 , ...,X n ,t) = u(xi,X 2 , ...,X n ,0,t) 

is the solution to the n-dimensional biwave equation for all \t\ < T. Since T is 
arbitrary, so u is the solution to the Cauchy problem in even dimensional space W 1 . 

Lemma 2.4. Given a function f : W l+1 — > K, which does not depend on the last 
variable, i.e. f(xi,X2,-~,x n +i)=g{xi,X2,-~,Xn), then 



M t (f)(x,0) = 



g(x + tz) 



dz. 



jb„(o,i) yl _ \ z \ 2 
Proof. Observe that, for x = (x, 0) and y — (y, y n +i), we have 



M t (f)(x) 



1 



^n+l JdB n+1 (0,l) 

We use the spherical coordinates given by 



f(x + ty)da(y). 



yi = smipi smip2— sin<£>„_ 2 sm<^„_i am tp n , 
y 2 = sin^i simp 2 ... sin<^„_ 2 sin<^„_i cos^„, 
y 3 = sin^i sin^2---sin^„_2COS^„_iCOS^ n , 



y n = sm^i cos ip 2 , 
y n+ i = cos 931, 

where < ipk < 7r, k = 1,2, ...,n — 1 and < <p n < 2ir. The Jacobian of this 
transformation is calculated as 



• n — 1 • n— 2 

sin if i sin ip 2 ...sm.ip n -i. 



Therefore 



M t (f)(x,0) 



1 



w n+l Jo 



g(x + ty)Jdcpidip2"-dip r , 
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Let us give r = sin ipi , and 



z\ = rsmtp2— sin</3„_x sin^m 
z 2 = rsmip 2 ... sin (p n -i cosip n , 

z n = r cos ip 2 . 

The Jacobian of above transformation is calculated by the formula 
J' = 



1 



r n 1 sin" 2 tf2 sin" 3 ip^... sin ip n -i 



Finally, we obtain that 
M t (f)(x,0) 



IjJn+l Jo 

7T />27T /•! 



7r /.2w rl 



W n +1 Jo 

2 



Jo 
(7(2; + tz) 



g(x + tz) 
1 



COS </5i 

Jdrdip 2 ...dip n 



-Jdrdip 2 ...dip T , 



Jo 

7T p2TT pi 



w n+l Jo 

2 



g(x + tz) — = JJ dz\dz% dz n 



Jo Jo 
g(x + tz) 



dz. 



/b„(o,i) V 1 ~~ M 2 
So the lemma is proved. 

We use the notation Ait(f)(x) = 



f{x + tz) 



□ 



dz, which is called 



^ra+l Jg„(0 1) a/1 - |z| 2 

as modified spherical mean of / (see e.g. [SS971 ISS03] ). Applying the result of 
Lemma 2.4, we obtain the formula of the solution to the biwave equation in the 
even dimensional space R": 

Theorem 2.5. The Cauchy initial value problem for the homogeneous n-dimensional 
biwave equation, where n > 2 is an even number, has the following solution 



u(x, t) 



1 



(n - l)!!(a 2 - b 2 ) 
b 2 d ( 1 a 



75R5 (W-'AMMM 



6 l&^di 



(ot) n - x A< t(^i)(a;) 



+ 



a /, 



bt 

t pau 



™) ((*)" )(,■]),/> 



Jb. 



s 9s 



((«) n - x A<.(0a)(a:) 



dsdv 



By a similar idea with the Duhamel principle for wave equations, the solution of 
the Cauchy problem for the nonhomogeneous biwave equation will be given at the 
next theorem 
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Theorem 2.6. The solution of the equation Ji.i|) - j"Op takes the form u = u + v, 
where u is the solution of the equation 12. l\) - ![27^) and 



v(x,t)— J ui (x,t,T)dr, 



where lj {x, t, r) is the solution of the homogeneous biwave equation 
d 2 oA fd 2 u 



dt 2 

with the initial conditions 
dui 
~dt 



a 2 A 



= 0, 



dt 2 
d 2 u 



IrA^j )=0, t>r, 
8 3 uj 



= 0, 



dt 3 



Proof. We start with the observation that 

^~(a 2 + b 2 ) ^Av + a 2 b 2 A 2 v 
t 

f(x,t) + 



I (^~( a2+b2 )w 2Auj+a2b2A2uj ) dT 



= f(x,t). 
Then, the above identity follows that 



Moreover, 



^-(a 2 + b 2 )^Au + a 2 b 2 A 2 u = f(x,t). 



0i (or) + = 0i (a;) , 



u \t=o -- 


= u \t 


=0 + z 




= 4> 


du 




du 




dv 






~dt 


t=o 


~di 


+ 

t=0 


di 


t- 


=0 


d 2 u 




d 2 u 




d 2 v 




dt 2 


t-Q 


~dt? 


+ 

t=0 




t=0 


d 3 u 




d 3 u 




d 3 v 




'dt 3 


t=0 


In 3 


+ 

t=0 


dt 3 


t=0 



02 (X) + = 2 (x) , 



03 (x) + = 03 O) . 



So the theorem is proved. 



3. Example 

Let us give an example demonstrating Theorem 12. II 



Consider the equation 



df 2 ~ dx~ 2 M dt 2 4 dx 2 



with the initial conditions 

du 



t=o 



dt 2 



COS X, 



t=0 



d 3 u 
dt 3 



□ 



(3.1) 



(3.2) 



t=o 
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The solution of the equation (|3.1[) - (|3.2j) is given by the formula 

«(*,*) = I (4cos (|) cos,-4cos;cos,- (-8sin(0 + sint) sin, 
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